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The Moody-Shapere-Wilczek's adiabatic effective ffamiltonian and Lagrangian method is devel- 
oped further into the matrix effective Hamiltonian (MEH) and Lagrangian (MEL) approach to a 
parameter-dependent quantum system. The matrix operator approach formulated in the product 
integral (PI) provides not only a method to find wave function efficiently in the MEH approach 
^— ■>) . but also higher order corrections to the effective action systematically in the MEL approach, a la 

' the Magnus expansion and the Kubo's cumulant expansion. A coupled quantum system of a light 

particle of harmonic oscillator is worked out, and as a by-product a new kind of gauge potential 
(Berry's connection) is found even for nondegenerate case (real eigenfunctions) . Moreover, in the 

QPI formulation the holonomy of the induced gauge potential is related to the Schlesinger's exact 
formula for the gauge field tensor. A superadiabatic expansion is also constructed and a general- 
ly^ , ized Dykhne formula, depending on the contour integrals of homotopy class of complex degenerate 
points, is rephrased in the PI formulation. 
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d . I. INTRODUCTION 

There has been much interest in a quantum system depending on a set ofparameters (parameter-dependent system) 
since the geometric phase for cychc evolution was discovered by Berry [ij and interpreted as a holonomy in a line 
bundle Q or a vector bundle [1, 0] over the parameter space. Since then, various approaches to the geometric 
^ ' phase have been introduced Q . Recently there have appeared some methods to find the quantum evolution operator 
(parameter-dependent vifave function) in the superadiabatic basis Q, improving the accuracy of the adiabatic basis 
■ and taking into account the transitions among the eigenstates of the Hamiltonian. An additional geometric factor 
, responsible for the transition amplitudes, known as the Dykhne formula 7], was found for a Hamiltonian system 
' having complex degenerate points @-[ll|. Joye, Kunz and Pfister have found a generalization of the Dykhne formula 
I I for a matrix Hamiltonian TT] . 
. The matrix effective Hamiltonian (MEH) and the matrix effective Lagrangian (MEL) were initiated by Moody, 
Shapere and Wilczek 9]. In the Hamiltonian approach, the quantum evolution operator, expanded in the eigenstates 
of the Hamiltonian operator, becomes a matrix differential equation 0, [l^. In the Lagrangian approach, the full 
time-evolution kernel, expanded in the basis of the Hamiltonian eigenstates, becomes an ordered matrix integration 
, 0. The MEL could be obtained by performing the ordered integration, but it has never been done due to the 
Jh ' noncommutativity of the matrix- valued operator on the parameter space; their formalism has not been completed due 
. . .' to the lack of a method to handle the matrix- valued operators systematically. 

On the other hand, independently of the MEH and MEL approach to a parameter-dependent quantum system, 
a method formulated in the product integral (PI) [13^ was introduced for the Wheeler-DeWitt equation, a time- 
dependent Klein-Gordon like equation, frequently encountered in quantum minisuperspace cosmological models, whose 
solutions were found explicitly using the properties of the PI [l^ . The positive and negative frequency decomposition 
was done by the method which can be regarded as a superadiabatic expansion method applied to relativistic quantum 
systems |15| . For a time-dependent harmonic oscillator, the quantum evolution operator and the corresponding 
Feynman kernel were found formally in the PI and expanded in ordinary integrals a la the Magnus expansion ^6] in 
Ref. [13. 
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In this paper we shall develop further the Moody-Shapere-Wilczek's adiabatic effective method into the MEH and 
MEL approach to a parameter-dependent quantum system. The PI formulation provides us a powerful method to 
treat the matrix-valued operators and also a unified picture to the MEH, MEL, geometric phases, superadiabatic 
expansion and geometric transition amplitudes. The matrix operator approach formulated in the PI provides not only 
a method to find wave functions efficiently in the MEH approach but also higher order corrections to the effective 
action systematically in the MEL approach, a la the Magnus expansion and the Kubo's cumulant expansion. The 
two-component wave function and the effective action are worked out for a quantum system of a heavy particle coupled 
to a light particle of a harmonic oscillator. In the MEL approach a new kind of gauge potential (Berry's connection) 
is found even for nondegenerate eigenstates (real eigenfunctions). As a by-product of the PI formulation it is observed 
that the holonomy of a nonabelian gauge potential (Berry's connection) is related to the Schleginger's formula for 
a gauge field tensor [l^, which is known as the nonabelian Stokes theorem [l9|. The Schleginger's formula is an 
exact formula for the geometric phase. A superadiabatic expansion is also constructed and a generalized Dykhyne 
formula, depending on the contour integrals of homotopy class of complex degenerate points, is rephrased in the PI 
formulation. 

The organization of this paper is as follows. In Sec. II the MEH for a parameter-dependent quantum system 
will be found, and a gauge transformation for the wave function and a unitary transformation for the Hamiltonian 
operator will be introduced that can separate the matrix-valued vector gauge potential from the Laplacian operator 
into the diabatic term. The MEH will be applied to a system of a heavy particle coupled to a light particle of a 
harmonic oscillator and the two-component wave function will be found in the PI formulation. In Sec. HI the MEL 
will be found using the spectral resolution for quantum system consisted of subsystems interacting each other. The 
matrix effective action will be constructed using the Magnus expansion and the Kubo's cumulant expansion. The full 
time-evolution operator kernel will be defined by a path integral using the matrix effective action. The MEL will be 
applied again to a heavy particle coupled to a light particle of the harmonic oscillator. In Sec. IV the geometric phase 
for a parameter-dependent Hamiltonian will be reformulated as a holonomy of a matrix- valued gauge potential. The 
higher order corrections to the Berry phase will be considered a la the Magnus expansion. The holonomy of the gauge 
potential will be related to the Schlesinger's exact formula for the gauge field tensor. In Sec. V a method for the 
superadiabatic expansion will be introduced in the PI formulation. In Sec. VI the geometric transition amplitudes 
for a quantum system having degenerate points coming from level crossing of eigenvalues will be shown to depend on 
the contour integrals of homotopy class around the degenerate points. 



II. MATRIX EFFECTIVE HAMILTONIAN APPROACH 
A. Formulation 

Let us consider a parameter-dependent quantum system whose time-independent Schrodinger equation is given by 

[Hi{R) + Hii{f,R)]^ir,R)=E^{f,R), (1) 

where the second Hamiltonian Hnif, R) depends continuously on i? as a parameter. This kind of quantum systems 
occurs frequently in many-body quantum system consisted of subsystems interacting each other. If ^{f,R) describes 
an atomic wave function for which Hi{R) is a sub-Hamiltonian for a nucleus (heavy particle) and Hii{f, R) another 
sub-Hamiltonian for an electron (light particle) , then the variable f is called the fast variable and the variable R 
the slow variable. The method we shall use below depends on the assumption that there exists a complete basis 
of the eigenstates for Hu (r, R) , allowing some possible degeneracies corresponding to level crossings of .R-dependent 
eigenvalues. This assumption is not too restrictive since the conditions for spectral resolution can be generally satisfied 
in most quantum systems. 

The Hamiltonian Hjj{r, R) will be assumed to have a complete basis of orthonormal eigenstates \ja{r, R)), 

Hiiir, i?)|ja(r, R)) = e,{R)\ja{r, R)) (2) 

with the orthonormality 

{kl3{r,R)\ja{r,R))^5,kSc.p, (3) 

where j, k run for continuous as well as discrete indices and a, /3 account for possible degeneracies. In this complete 
basis the wave function 4'(r, R) can be expanded by 

*(r, R)^Y. = *^(^: m{R). (4) 

ja 
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where we have introduced a compact notation, arranging the basis \ja{f,R) as a column vector $(r, i?) and the 
ampHtudes ipja{R) as another column vector ^'(i?), and denoted the transpose of matrix or vector by T, and assumed 
the usual matrix operations. 

Thus, the Schrodinger equation may take the following matrix form 

[<^*{r,R)Hi{R)<^'^ {f,R) + Hii,d{R)]-^{R) = EI^{R), (5) 

where * denotes the dual (bra) vectors and Hjj jj(R) is a diagonal matrix of eigenvalues ej(R) arranged in the definite 
order used for the column vectors $(r, i?) and and / is the identity matrix. The matrix operator in the left 

hand side is a MEH. For a Hamiltonian Hj{R) of the form 

Hj{R)^^ + Vi{Rl (6) 

we introduce a matrix- valued gauge potential [23| by 

A{R) - i$*(r,i?)V^<J>^(f,^), (7) 

in components Aja^kfsiR) — *Ock('^i -R)|V^|fc/3(r, i?)}. Then the Schrodinger equation (in unit of ?i = 1) separates 
completely as 

-^(IWj^-tA{R)f +IVi(R)+HiiMR)]'^iR)=EI^{R). (8) 

The term Hjjj:){R) represents a back- reaction of the system Hii{f, R) on the system Hj{R). In the Born-Oppenheimer 
approximation in which one sets A{R) — 0, Hu oiR) contributes to the effective potential 

V.siR) = IVi{R) + Hi I MR)- (9) 

However, this is not the only effect that the system Hjj{r, R) has on the system Hj{R) because it also introduces the 
matrix-valued gauge potential to the kinetic energy term. The induced gauge potential has an important quantum 
mechanical effect, that is, the Berry phase Equation ([8]) has been written as a compact form for a vector 

equation that was already derived in Refs. [20, l2l|. What has been obtained hitherto is the fact that the original 
parameter-dependent Schrodinger equation separates as a vector equation and the Hamiltonian operator becomes a 
matrix operator called MEH operator 

H.siR) = ~^{I^R - ^AR))^ + IVi{R) + HuAR)- (10) 

More generally, the MEH operator takes the form 

H,s{R) = ^*{r,R)Hi{R)<^^{r,R) + HuAR). (H) 

where the first term is a matrix of expectation values with respect to the basis \ioi{r, R)). 

We may separate the gauge potential from the Laplacian operator in the parameter i?-space by first expanding the 
wave function in the complete basis and then inserting a gauge factor defined in terms of the PI: 

R 

^{f,R) = <S>^{f,R)(j[exp[iA*{R) ■ dR\^^{R), (12) 

Here we have introduced a PI for an operator- valued function (the matrix- valued gauge potential ([7])) to replace the 
ordered integral which is a symbolic notation in quantum field theory; however, in our opinion the PI should deserve 
much more attention than the ordered integrals in that the former has a more rigorous mathematical foundation (T3| 
than the latter. Finally, the Schrodinger equation becomes 



J R _^ R 

^—\/l + IVi{R) + (jlexp[iA*{R)-dR\) HiiAR){Y[eM^A_*{R)-dR\)]^iR) = Eli'{R). (13) 



The gauge potential A{R) in Eq. (fTO|) has been eliminated and the matrix Laplacian operator has become a diagonal 
matrix Laplacian operator in Eq. p^ . whereas the back-reaction of the system Hij{f,R) on the system Hj{R) has 
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gained the off-diagonal matrix (the last term), which is known as a diabatic term in literature [22|. It should be 
noticed that Eq. (jl2p is a gauge transformation of the original wave function and Eq. (jl3p is a unitary transformation 
for the MEH. There appears a similar singular gauge transformation in anyon system with the Chern-Simons action 
[23} : however, the similarity is expected naturally because the singular gauge transformation of anyon system is a 
kind of the gauge transformation (IT2]) with an abelian singular gauge potential A{R) and the PI becomes an ordinary 
integral: 

Jl exp [iA* (R) ■ dR\ = exp [i / A* {R) ■ dR] (14) 

In the anyon system the boundary condition for the original single- valued wave function ^(r, R) is transferred into a 
different boundary condition for the gauge transformed wave function vE'(i?), for example, such as the multi-valuedness. 

The ultimate goal of the MEH approach to the quantum system ([T]) is to find the energy eigenvalues and the wave 
functions. To find the energy eigenvalues, we need to solve the MEH equation (jl3p . which involves an unknown energy 
eigenvalues E. A typical method would expand the vector-valued wave function ^{R) by some basis of orthonormal 
functions of R. The MEH equation becomes a blocked matrix eigenvalue problem. The wave function can be obtained 
as an eigenvector; however, there is a more efficient method that rewrites the MEH equation (|13p by a two-component 
equation of the first order. This method was successfully applied to find the complete set of the wave functions for the 
quantum Friedmann-Robertson- Walker cosmology minimally coupled to scalar fields in Refs. [l^. [l5|. The explicit 
exhibition of this method to determine the wave functions will be done for an example. 



B. Applications 

As an application of the MEH approach to quantum evolution operator, we shall consider the Hamiltonian for a 
system of a heavy particle coupled to a light particle with the potential of a specific form 

Hir, R)^^ + Vj{R) + ^ + ^u^'iR^. (15) 
ZM Zm Z 

Equation (jl5p is an effective Hamiltonian for two interacting particles, because if M is the mass of the heavy particle 
and m the mass of the light particle, then the potential terms describe a small perturbation about the instantaneous 
position R. The super-Hamiltonian of the Wheeler-DeWitt equation for the Friedmann-Robertson- Walker cosmology 
minimally coupled to a free massive scalar field has the form (jl5p . in which R corresponds to the scale factor of 
the universe while r corresponds to the scalar field [13, [Hi • The light particle of the Hamiltonian ((T5|) describes a 
three-dimensional harmonic oscillator with a variable frequency depending on the position of the heavy particle. The 
exact wave function and energy eigenvalue for the Hamiltonian (|15l) can be determined by a limiting process in which 
the light particle acts on the heavy particle as the back-reaction and the heavy particle acts on the light particle 
through the variable frequency, and so on. 

To simplify the model further without any loss of generality, we shall restrict to the two-dimensional problem of 
the heavy particle coupled to the light particle of a harmonic oscillator 

Hi.,X) = — + V,iX) + l- + -.\X).^. (16) 

The orthonormal eigenfunctions of Eq. ([2]) are the harmonic wave functions 

= {^0^y^'H,{V^^^)x)e^p[-lmu;iX)x'), (17) 

where Hj{x) are the Hermite polynomials. It is convenient to work with the number-state representation and the 
creation and the annihilation operators act on it. The creation and the annihilation operators obviously depend on 
X as a parameter 



a\X) 



mujiX) 1 

X ^ i — , -P, 

2 y^2moj{X) 



, > mujiX) 1 
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The number states depend also on X as a parameter and the harmonic wave functions ([T7|) are the coordinate 
representation of the number states. It is straightforward to see that the Hamiltonian of the hght particle in the 
number-state representation leads to the diagonal matrix 



HiiMX) ^ i^{X)(a^ {X)a{X) + 
and so does the matrix- valued gauge potential ([7]) 

duj{X)/dX p 2.^x ^2fY^l 



(19) 



(20) 



It should be remarked that for the real eigenfunctions ([TTl) the matrix- valued gauge potential ([20| does not contain 
any diagonal part. The operators occurring in Eqs. (fTQ]) and (|20|) 



Ni = a''{X)a{X) + 



I 



a^'{X) + a^{X), iVa = aTi(x) -a^(X), 



(21) 



constitute a Lie algebra so(2, 1) with the basis {iVi, A^2/2, A^3/2} satisfying the commutators 



N2 Ns 
2 ' 2 



27Vi, 



(22) 



In the instantaneous number-state representation the operators pip represent constant matrices, and the constancy 
of these operators will be understood throughout this paper. 

With the aid of commutators ([22|) the diabatic (last) term in Eq. (fTS]) turns out 



X 



X 



II. D 



Y[exp[iA*{X) -dX] 

ldu:{X)/dX^ 



HH.D{X)\J{e^I>[iA*{X) ■ dX] 
1 . . (d^{X)/dX^ 

L0{X) ) 



Then the MEH equation becomes 



'2M dX^ 



IVi{X) + Hii,d{X) §(X) = EI-^{X). 



(23) 



(24) 



The diabatic term for the back- reaction of the light particle consists of both the diagonal entries (the first term) and the 
off-diagonal entries (the second term) in Eq. ((23|) . Therefore, Eq. p4| is not the standard energy-eigenvalue problem. 
We may expand the column vector wave function '^{X) by some basis of orthonormal functions of X; for example, by 
the eigenfunctions for the heavy particle, Hj{X) = P'^ /2M + Vi{X). Then Eq. turns out to be a blocked matrix 
energy-eigenvalue problem, whose entries are again matrices of the same dimension as the eigenfunctions for Hj(X), 
which may be infinite dimensional. This causes a great difficulty in finding the energy eigenvalues. We may circumvent 
the infinite dimensionality problem by truncating the blocked matrices to an appropriate size or by rearranging the 
double indices, for instance, according to the sum of both indices. 

Let us turn to the problem to find the wave function for the Hamiltonian (|16l) . We can make use of the eigenfunc- 
tions determined from the matrix equation ([24]); in fact, this procedure involves the double expansion by harmonic 
wave functions of x and the eigenfunctions for Eq. (|24|). Thus the resulting wave functions may serve as a good 
approximation in some regions of the configuration space of x and X. These wave functions do not, however, show 
the asymptotic forms. Instead, we may employ the two-component wave functions. The detailed calculation is shown 
in Appendix [BJ The final result is the two-component wave function 



'^{x,X] 

d^{x,xj 
dX 



Xo 







<p'^{x,X) 



iA*{X) I 
2M{IVi{X) + Hii^d{X)- EI) iA*{X) 



dX 



7$(a;',Xo)*(x',Xo)dx' 
^'jci>(x',Xo)^^5^dx' 



(25) 



What the two-component wave function implies is that once the energy eigenvalue is determined by some methods, 
the wave function at an arbitrary position of the heavy particle is evolved from the wave function at a given position of 
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the heavy particle according to Eq. ([25]) . Of course, the two-component wave function for the Hamiltonian ([16]) does 
not belong to the Cauchy initial value problem but belongs to the elliptic-type eigenvalue problem. Nevertheless, there 
is still some similarity in finding the wave function between the Hamiltonian ([T6| and the super-Hamiltonian for the 
Friedmann-Robertson- Walker universe minimally coupled to a free massive scalar field as mentioned earlier; the only 
difference is that the latter has the opposite sign for the kinetic energy operator, indicating that the Wheeler-DeWitt 
equation is of the hyperbolic-type and the two-component wave function is truly the evolution of the Cauchy initial 
data. 

To know the wave function explicitly, one should be able to evaluate the PI in Eq. (|25p in terms of the known 
algebraic operations. There is no general method to evaluate the PI in a closed form except for a few cases. So we 
may use the useful properties such as (|A8|) and (|A10I) of the PI to extract the dominant terms and to expand the 



remaining part in a perturbation series. Some methods to evaluate the PI have been proposed in Refs. [IJ, [15j and 
used to find the wave function for the Friedmann-Robertson- Walker universe minimally coupled to scalar fields with 
various potentials. According to the result of Ref. pL^, the wave function ((25)) should take the form 

^{x,X) = X) [Pi {X)Qn{X, Xo) + P2 {X)Q2i{X, Xo)] Wr-^ [Pi (Xq) , P2 (Xq)] 



dXo 



Hx',Xo)-i'{x',Xo)dx' -P2{Xo) 



'^{x ,Xq} — dx 

CAq 



+$^(a:,X)[Pi(X)Qi2(X,Xo 
dPiiXo 



P2iX)Q22iX,Xo)]WT-^[PliXo),P2iXo)] 



J ^{x',X^)^{x',X^)dx' + P,{X^) J $(x',Xo)^^^-^dx' 
where Pi{X)^s are the fundamental solutions to the diagonal matrix equation 

+ IVi{X) + HhMX) - El] Pi,2(X) = 0, 
and Wr-i[Pi(Xo),P2(Xo)] is the inverse of the Wronskian matrix, and 

X 



Qll{X, Xq) Qi2{X,Xq 
Q2l{X,Xo) Q22{X,Xo, 



J|exp 



Lii(X) ii2(X)\ 

L2l{X) L22{X) 



dX 



where 



Lii{X) 
Li2{X) 



dX 

dP2{X) 



dX 



m-\Pi{x),P2{x) 

Wr-i[Pi(X),P2(X) 
i2i(^) = Wt-\Pi{X),P2{X)]\Pi{X)iA*{X) 
L22{X) = Wi-\Pi{X),P2{X)]\Pi{X)iA*{X) 



'^^'^^\a*{X)P,{X) P2{X)^A*{X)'^^'^^^ 



iA*{X)P2{X) - P2{X)iA*{X) 



dX 
dP2{X) 



dX 



dPi{X) dPi{X) . 



dX 
dP2{X) 
dX 



dX 



iA*{X)Pi{X) 



dPi{X). 



dX 



iA*{X)P2{X) 



(26) 



(27) 



(28) 



(29) 



The matrices Qij(X, Xo)'s can be determined by using the Magnus expansion (jA5| or the T-ordered integral. This 
(superadiabatic) expansion of the PI was used to find the asymptotic forms of the wave functions ^15*| . 



III. MATRIX EFFECTIVE LAGRANGIAN APPROACH 



A. Formulation 



In the previous MEH approach to quantum evolution operator, one should solve the vector equation (I13p . which is 
equivalent to the parameter-dependent quantum system ([1} through the spectral resolution and the gauge transfor- 
mation (fT2|) . In order to use the MEH approach, however, one should surmount the difficult task to find the exact 
solutions to Eq. (|13l) . In this section we shall develop the MEL approach, the second half of the matrix operator 
approach, to quantum evolution operator. 

Let us consider again the following Lagrangian for the quantum system ([1} 



L{r,R)^Li{R)+Lii{r,R) 



(30) 
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Due to the parameter J?-dependence of the Lagrangian Ljj{r,R), a common method of path integral performs the 
r-path integral first, whose result is an i?-dependent operator; on the other hand, the -R-path integral might be very 
difficult to perform in general without perturbation techniques applied to the i?-dependent operator. Contrary to the 
path integral, in the MEL approach we shall apply the spectral resolution of the system // to the full time-evolution 
operator (propagator in the configuration space (r, R)) and perform the remnant matrix operator path integral for the 
parameter _R-space. The MEL is a nonperturbative one giving the i?-dependent operator a matrix operator, which is 
not necessarily diagonal. 

The full time-evolution operator kernel can be written as a Feynman sum over all the paths from a configuration 
{fi,Ri) at time ti to another configuration {rj^Rf) at time tf 



K{ff,Rf,tf;fi,Ri,ti 



^1 



D[R] D[r\e^p i (Li{R) + Lii(f,R))dt 



(31) 



Noticing that the total Lagrangian operator depends on r through Lu^r, R) only, we may first perform r-path integral. 
The quantum state of the system // in the basis of eigenstates of the Hamiltonian Hjj (r, R) propagates as 



Un{ff,Rf,tf;'ri,Ri,ti) = 



J D[r\exp^i J Lii{f,R)dt 

'^^irf,Rf)Ylcx^p[~i(^Hii^DiR) +^iR)- ^*{f^,R^) 



(32) 



The parameter R depends on t implicitly so that R — dR/dt. Then the full time-evolution operator kernel in the 
spectral resolution of Hii{r, R) takes the form 

K{ff,Rf,tf;n,R,,U) = <i>'^{ff,Rf) D[R]cyi^\i j\i{R)dt 

*/ 

X exp [-1 [Hu,d [R) +^{R)- R) di] ^*{r^,R^). (33) 

ti 

Next, we can factor out the diagonal matrix from the PI using the general perturbation formula (sum rule) (Tsj 



W exp \~i {Hii,d {R)+A^{R)-R 



dt 



nexpUe*"(^'^''A5(i?)-7?e^ 



in(R,Ri 



dt 



(34) 



where we have defined a frequency matrix by 



(Hn,D{R) + ^D{R)-R)di 



(35) 



where A£,{R) and Ao{R) denote the diagonal and the off-diagonal parts of the matrix- valued gauge potential. This 
separation of the frequency matrix makes the Magnus expansion of the ordered integrals converge in the adiabatic 
regions, avoiding the possible divergence noticed in Ref. |24| . 

In the MEL approach, we may write the full time-evolution operator kernel p3p involving the PI by an ordinary 
exponential operator a la the Magnus expansion (see Appendix |X1 especially Eq. (lASp ) 



K{ff,Rf,tf,n,R„U) = <S>'{ff,Rf) D[R]e^p i I {Li{R) ^ Hn^oiR) ~ AI{R) ■ R)dt 



X exp 



M[-te'"(^'^^^I^{R) ■ Re-'"M)dt\ ^*{r,, R,). 



(36) 



The two matrix exponentials do not commute each other; therefore, we need the Baker-Campbell-Hausdorff (BCH) 
formula to get a single matrix exponential. But the powerful cumulant expansion method by Kubo [25] extracts some 
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pieces of information about the evolution operator without the BCH formula. To obtain the diagonal matrix element 
of the evolution operator, it suffices to find the diagonal matrix element of the second matrix exponential 



exp 



exp 


(exp 











(^57) 



where the average (• • • ) is defined to take the diagonal element, and the subscript c denotes the cumulant expansion of 
the matrix operator (see Appendix [BJ . Then, the matrix effective action results from the diagonal part of exponential 
operators in Eq. (|36|) 



{Li{R) - HuAR) - Al{R) ■ R)dt + (exp[ / M[-ie'^^^-'^^^ Al{R) ■ Re-'^^^^^^^Ut 



(38) 



The full time-evolution operator kernel up to the diagonal matrix effective action is approximated by 



K{rf,Rf,tf,n,R„U)^<^' {rf,Rf) / i?[i?]exp iS^s[R) '^*{n,R^)■ 



(39) 



It should be pointed out that the time-evolution operator kernel using the matrix effective action is not exact since 
the cumulant expansion of the matrix operator was undertaken with the diagonal matrix average in Eq. p7p . To 
include the transition effects it needs to insert a transition matrix of the form I + T, T being an off-diagonal matrix 
that was not included in the diagonal matrix average. 



B. Application 



As an application of the MEL approach to quantum evolution operator, we shall consider again the system ([IB] 
whose Lagrangian is 



L{x, X) ^ fx' - Vj{X) + - y-^(X).^ 
In order to find the full time-evolution operator kernel (j36l) . we need to perform the Magnus expansion 



(40) 



exp 



(41) 



With the choice of the real eigenfunctions p7|) for the light particle, the matrix-valued gauge potential (PO]) docs not 
contain a diagonal part so that 



m.X.)^ j\{X)dtN,^ A^(X) = ^^^;^^A^3. 
Using the Lie algebra and the BCH formula, we can show that 



(42) 



exp 



exp 



M[X{h{X)N:, + if2{X)N2)]dt 



(43) 



where 



/3(^) 



dLu{X)/dX 
4w(X) 



cos(2 f uj{X)dt^, f2{X) 



duj{X)/dX . 



4a; (X) 



sm(2 / Lu{X)dty 



(44) 
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The first few terms of the Magnus expansion (jA7l) are 
Af(o) = X{h{X)N3 + if2{X)N2), 



where 



X 



M(i) = -4iX / dXiF[,o]{Xi,X)Ni, 

Af(2) = ^iX [ dXi [ ' dX2F^^^{X2,Xi){MX)N2 + ih{X)N^) 



Xi 



Xi 



X 



M, 



(3) 



Xi 

-16iX [ dXi 
Jx, 



X 



-iX 1^ dXi dX2F[,o]{Xi,X){MX2)N2+if2{X2)N3), 



Xi 



dXo 



Xi 



X2 



Xi 



dX3F[32]iX3,X2)if3iXi)hiX) + h{Xl)f2{X))Ni 



16 . • 



X 



Xi 



Xi 



Xi 



+ — lX I dX, I dX2 I dX3F[2i](X2,Xi)(/3(X3)/3(X) + /2(X3)/2(X))iVi, 



Xi 



= f2{x,)hix,) - Mx,)f2ix,). 

From the Lie structure the Magnus expansion must be closed in the form 

rtj 



exp 



exp 



where G^'s are determined by summing over M(i)'s, that is, 

XJ2G^{X)N,^Y.^h^) 



dtX{Gi{X)Ni + G2{X)N2 + Gs.{X)N3) , (47) 



(45) 



(46) 



(48) 



Therefore, the fuh time-evolution operator kernel (1361) becomes 



Xf 



K{xf,Xf,tf-x,,X,,U) = ^'{xf.Xf) I D[X]e^p i / (LiiX) - uj{X)Ni)dt 



Xi 



X exp 



dtX{Gi{X)Ni + G2{X)N2 + G3{X)N3) <^*{x,,X,). 



(49) 



Here $(x,X) represents the column vector of the eigenfunctions (|17p of the harmonic oscillator. 

The matrix effective action is obtained from Eq. (1551) by evaluating the Kubo cumulant expansion (jCl[) in which 
the average (• • • ) is taken for the diagonal element. After some algebra, we find the first few terms of the cumulant 
expansion 



exp 



dtX{Gi{X)Ni + G2{X)N2 + G3{X)N3) 



dtXGi{X)Ni 



1 



dtXG2{X)\ - / dtXG^iX) 



dtXGiiX))( / dtXG2 



j(X))'7V4 + ' dtXGiiX)'^ ' dtXG3(X)) + • • • , (50) 



where iV^'s are the constant matrices in the instantaneous number-state representation for the operators 
iV4 = a^^{X)a^X) + a^X)a^^{X), N5 = a^HX)a^{X) - a^X)a^\X). 



(51) 



Finally, up to the lowest Magnus expansion and the Kubo's cumulant expansion we obtain the matrix effective action 
for the heavy particle 



S, 



off 



dt 



(52) 
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where 



Introducing the mechanical momentum P = MX for the heavy particle, the matrix effective action ([5^ becomes 
S,s = "^^i^i^ - 9{X)N,f - Vi{X) - cj{X)m - ^9HX)N!] . (54) 

One can easily recognize that the effective action for the heavy particle involves a gauge potential g{X)Ni and 
an additional back-reaction g'^{X)Nf /2M besides the back-reaction uj{X)Ni from the light particle. The gauge 
potential results from a sequence of transitions accompanying the creation of two quanta followed by the annihilation 
of two quanta or vice versa in the instantaneous number-state representation. Contrary to the common belief that 
the induced gauge potential does not exist for the nondegenerate case (real eigenfunctions) , the MEL approach to 
quantum evolution operator shows that up to the lowest Magnus expansion and the Kubo's cumulant expansion 
the matrix effective action has manifestly a new kind of mode-dependent gauge potential (Berry's connection). For 
the lower modes the additional back-reaction g'^{X)Nf /2M is suppressed by m/M compared to the back-reaction 
a;(X)A^i, but for the higher modes it is enhanced by Nim/M. 



IV. GEOMETRIC PHASE AND SCHLESINGER'S FORMULA 

Let us consider a Hamiltonian H{R) depending continuously on a parameter-vector R in the parameter space which 
needs not to be of the specific form ([T]), but just depends on i? as a parameter. In the spectral resolution, for each 
fixed R the Hamiltonian satisfies the eigenvalue equation 

F(^)|ja(i?))-e,(i?)|ja(i?)), (55) 

where the Latin letters j, k, - ■ ■ denote the energy levels and the Greek letters a, f3, - ■ ■ denote the possible degeneracies 
of the Hamiltonian. As the Hamiltonian H{R) changes continuously along a curve C{R), the basis of the eigenstates 
\ja{R)) changes also continuously along the curve, which has an analog of a moving coordinate system in classical 
mechanics. We may arrange the eigenstates \ja{R)) as a column vector <I>(i?) as in Sec. IhI and introduce again the 
matrix- valued gauge potential [1, US] ^{R) (or the coupling matrix [17]) from Eq. ([7]) in the form 

V^$(i?) = ~iA{R)'S>{R). (56) 

Equation ([56l) is one of the Cartan's structure equations for a vector bundle over the /?-space, which is the torsionless 
condition relating the change of the orthonormal basis $(/?) to the matrix- valued connection A{R) [2^. Under a 
change of frame, ^{R) = U{R)^{R), corresponding to another choice of basis, it can be shown that the gauge potential 
transforms truly as a gauge field 

A{R) = U{R)A{R)U''{R)+i{Vj^U{R))U\R), (57) 

rather than as a tensor. It is straightforward to see from the orthonormality of the basis that the gauge potential 
is Hermitian, A^(i?) — A{R). The gauge potential may be divided into the geometric part A'^{R) coming from the 
arbitrariness in the choice of the phases in degenerate eigen-subspace of the Hamiltonian and the dynamical part 
A^{R) such that 

A{R)=A^{R) + A^{R), (58) 

where 

ej(i?) - ek{R) 



Here and hereafter the superscript D is used for the dynamical parts in this section while the subscript D was and will 
be used for the diagonal matrix parts throughout the paper, but this would not cause any confusion. The geometric 
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part A'~^ (i?) can not be determined from the eigenvalue problem ([55]) and is the nonabehan effective gauge potential 
[20| . also known as the Berry's connection m. 

The transport of the frame of eigenstates along a curve C{R) beginning at a base point Rq is given in the PI by 



C,Ro 



(60) 



The nonabelian gauge field tensor for a vector bundle of the gauge potential over the i?-space follows from a commu- 
tator acting on the orthonormal frame ([BH]) as 



d 



d 



dXra ' dX 

n 



.fdAcAR) dAc,m{R) 



dXn 



- [Ac,™(i?),Ac,„(i?)] ]^c{R) = F^^{R)^c{R), (61) 



where the derivative d/dXm in the i?-space is actually a covariant derivative in the horizontal subspace of the tangent 
space of vector bundle [26] . Note that the nonabelian gauge field tensor Fmn {R) for the induced gauge potential ([55)1 
satisfies the same field equation as a true gauge field, and so the field equation allows the interpretation of d/dXm as 
a covariant derivative. For an abelian gauge potential, the commutator is the curl operator and can be used for the 
(abelian) ordinary Stokes theorem. Equation (I61|) relates the curl operator to a gauge field tensor for a nonabelian 
gauge potential too. The nonabelian Stokes theorem was discussed in Ref. [l^. However, in the history of the PI (or 
matricant) the nonabelian Stokes theorem dated back to the Schlesinger's formula [31 that related the holonomy of 
the affine connection to the surface integral of the Riemannian-Christoffel curvature tensor. The field tensor explains 
how the geometric phase for closed loops 0, 01 enters the evolution of parameter-dependent quantum system. 
The holonomy matrix in the PI ([501 may be evaluated by the Magnus expansion [16] as in Ref. 17 1 



n 

C(i?) 



-iAc{R) ■ dR 



= exp 



M[~iAc{R)] 

CiK) 



(62) 



whose first term only with the geometric part is nothing but the Berry phase. With the ordinary Stokes theorem 
applied we see that 



exp 



-I I A'g{R) ■ dR 
'c 



exp 



i / {VjixA2{R))-dS 
Js 



exp 



F§{R)-dS] 



(63) 



where dS is an area vector for the surface bounded by C{R) and F§{R) is the Berry's curvature (|6ip along that curve. 
Higher order corrections to the Berry phase may be calculated from Eq. (see Appendix |X] for higher order terms 
for the Magnus expansion.) 

We now turn to the Schlesinger's formula for a gauge field tensor. The holonomy matrices (p^ of the gauge potential 
for loops C{R) at the base Rq form a holonomy group at Rq. The Schlesinger's formula connects the holonomy of a 
gauge potential to a gauge field tensor. After the affine connection and the Riemannian-Christoffel curvature tensor 
of differential geometry are translated into the gauge potential and gauge field tensor of gauge theory, the holonomy 
(p^ is related to the gauge field tensor by a mixed product integral 



C{R) 



-tAc{R)-dR]= n cxp[£ fuF{^,v)F,p{^,r,)U^'{tv)j(^^)d^dr, 



(64) 



where (pp is a surface map from a rectangle {(x, y)|a < a; < 6, c < y < d} in a Cartesian coordinate system to the loop 
C{R) and its interior in the i?-space, j(^(f>p, <j)q/^,r]j the corresponding Jacobian, and 



Uf{£„v) ^Yi'^^P -'iAm{R{a,y)) 



, dfjjmia, y) , 
dy 



dy\ J|exp 



-zA„(i?(x,ry))^^^^dx 



(65) 



The Schlesinger's formula ([M)) is an exact one connecting the holonomy of a gauge potential to a gauge field tensor. 
In fact, the Magnus expansion for the holonomy ((62|) is one of the expansion method to the Schlesinger's formula. 

Using some of the above results obtained in Sees. HIl and Hill we shall develop the superadiabatic expansion in the 
next section, and for that purpose a time-dependent Hamiltonian with an adiabatic small parameter will be quite 
illustrative. 
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V. SUPERADIABATIC EXPANSION IN NONDEGENERATE CASE 



In this and next sections we shall reformulate the quantum evolution of a time-dependent system in the PI. Then 
the parameter is the explicit time and the parameter space is a real axis or a complex plane of time. The advantage 
of the PI formulation is some useful properties and the efficiency of the PI as a computational tool. 

In order to develop the superadiabatic expansion, a generalization of the adiabatic expansion, we shall insert 
explicitly an adiabatic parameter into the Schrodinger equation for a time-dependent Hamiltonian 



(66) 



where e is an adiabatic small parameter depending on the Planck constant h. The quantum evolution operator for 
Eq. (|66|). which evolves a quantum state at a time tq to the quantum state at an arbitrary time r, may be given in 
the PI by 



li^ir)) =l[exp\~-Hir)dr\\i;iTo)) 



(67) 



The superadiabatic expansion for the Schrodinger equation (|66p was worked out for a time-dependent harmonic 
oscillator in Ref. [l^l- In general, when expressed in the basis of H{t), the PI becomes 



Y[exp --H{T)dT ==$^(T)]^exp -ii-H{T)~A^{T)\dT $*(to). 



(68) 



where T denotes the transpose and * the dual operation. We can apply the general perturbation formula (sum rule) 
[l3{ to factor out the dynamical phases: 



]Jexp -i{-H{T) - A^{T)^dT 



exp 



Hi:){T)dT\ JJ^exp iAi{T)dT 



where 



A/(T)=exp - / HD{T)dT ^^(T)exp 



HD{T)dT 



Next, we shall use the formal equivalence between the PI and the time-ordered T-exponential (T^ - [l7| 

T CO .r .T-i .,.("-1) 

JJexp iyl/(T)dr =/ + ^r / Ai{T')dT' Ai{T")dT"--- Ai{T^"-^)dT^"-\ 

To n=l "'■^0 -Jto Jto 

and integrate by parts the rapidly oscillating factors in the adiabatic regimes to obtain 

r oo 

J|exp[iA/(T)dTj = / + ^e"A/,„(T,To). 



(69) 



(70) 



(71) 



(72) 



When the series in power of the smallness parameter e terminates at the iVth step, this gives another derivation (not 
exactly the same) of the Berry's superadiabatic basis 



\^n{t)) = $^(T)exp 



N 



HD{T)dT + ^""^lAr, To) **(To)|^(ro)) 



(73) 



It is worth to note that the superadiabatic expansion is well-defined and converges sufficiently in the adiabatic regimes 
in contrast to the Berry's diverging basis for large N. Also note that the superadiabatic basis (|73p resulted from the 
T-exponential (ITT]) , but this is not the only possible form for the PI. For example, the formal solution in the Magnus 
expansion (Appendix \^ may be employed to give the exponential correction 



IV'(t)) = ^'^{t) exp 



/ HD{T)dT exp^e"i,,„(r,To) $*(to)|V(to)). 



(74) 
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and the Nth truncation 

I^jv(t)) = $'^(T)exp 



N 

Hd (r ) dr] exp e"i7,„ (r, to )] $* (to ) I V (to ) ) . 



(75) 



Especially, when the group structure of the Hamiltonian operator is known, the Magnus expansion proved quite useful 
as in Ref. fTTj. Now, we turn to the case when there are complex degenerate points coming from level crossings. 



VI. GEOMETRIC TRANSITION AMPLITUDE 



Let us suppose that the Hamiltonian H{t) has an analytical continuation in a strip of a complex T-plane, denoted 
by z including the real T-axis, such that H{t) has a basis of eigenstates with discrete eigenvalues except at finite 
isolated complex points of the degeneracies. These conditions were used to derive the geometric transition amplitude 
for a two- level system and the applicability to general systems was suggested in Refs. 04111 ■ this section we shall 
reformulate the geometric transition amplitude formula in a general form using the analytic properties of the PI. 

If there are no complex degenerate points for the extended Hamiltonian H{z), then the path-independence of 
quantum evolution operator is a consequence of the identity operator for the PI for a closed loop 



n ^'^P 



-i(^^HDiz) - A^{z)]dz 



n 

Cll{z),To 



-i{--HDiz) - A' {z)]dz 



(76) 



for two different paths in the z-plane with the same initial and final points on the real axis. This path-independence 
makes the use of complex z-plane legitimate. However, the situation is not so simple when H{z) has degeneracies 
at isolated complex points since A'^{z) becomes a meromorphic operator, whose regular singular behavior can be 
seen from Eqs. ([55]) and ([5^ as the two or more energy eigenvalues approach each other. For the sake of simplicity 
we shall assume that a pair of complex conjugate points zq and ^o are the only degenerate points. Then any curve 
connecting the initial and final points is characterized by the homotopy class of winding numbers around zq or Zq or 
both Zq and Zq. We shall restrict the study to the homotopy class of the winding number around zq, but the choice of 
the homotopy classes should be determined from the preference of exponentially small transition amplitudes for the 
superadiabatic expansion. 

The quantum evolution operator for a curve C(z) consisted of Cq{z) and c["'\z) of winding number n, (C(z) = 
Co(z) + cf)(z)), is 



n *^^p 

C(z),n, 



~-H{z)dz 



n *5^p 

Co(z),To 



H{z)dz\ exp — H{z)dz 



(77) 



where H{z) — Hi:,{z) — eA^ {z). The Pi's for two different closed loops with the same winding number are equivalent 
up to a unitary transformation [l^ 



n 

C["\za) 



exp 



H{z)dz Yl e^P 



c['-"Hzo) 



--H{z)dz 
e 



(78) 



for some unitary matrix U . The unitary transformation equivalence of the ordered loop integrals for the winding 
number one was also discussed in Ref. [ll| . 

It is remarkable that the quantum evolution operator of the Hamiltonian with complex degenerate points has an 
additional factor determined entirely by the PI over the loops according to the homotopy class of degenerate points, 
which is geometric in nature. It should be noted that our derivation of the geometric factor was based only on the 
extension of the Hamiltonian to a strip of the complex time plane including the real time axis and the existence of 
isolated complex degenerate points. 

The PI of meromorphic operators over loops beyond the simple poles is technically difficult to perform, but a 
method may be put forward to replace the PI by the Magnus expansion (see Appendix E)) . Thus the PI over the loops 
becomes 



n 

c["\zo) 



exp 



'-H{z)dz 



exp 



M[~-H{z)]dz 



(79) 
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where the Magnus global expansion up to the third order is 

M[--H{z)] = -tH{z) + ± r d,'[H{z%H{z)] + -^^ f dz' f dz' [H{z'),[H{z H{z)]] 



' Zq J Zq 

The first term in the right hand side yields 



^ I dz' I dz [[H{z ),H{z%H{z)]] + --- . (80) 



exp 



I 



H{z)dz =exp -^^ — —H{zo) , (81) 



where H{zo) denotes the residue operator at zq. The geometric factor (|8T|) gives a generalization of the Dykhne 
formula for the transition amplitudes for the homotopy class. Equation (|79p includes higher order correction for 
the geometric transition amplitudes. A complex curve starting on the real axis and ending on the real axis can be 
homotopically deformed into a loop around one of the complex degenerate points and a segment of the real axis if the 
curve goes initially outside them. The evolution operator composes thereby of that of the real axis and a loop integral 
with the winding number one around that complex degenerate point. The loop integrals with the winding number 
one were discussed in literature [ll| . The PI approach in this section is more systematic and results thus obtained are 
more comprehensive than others, leaving the investigation of all homotopy classes and the choice of paths ensuring 
the convergence of the evolution operator for future research. 

VII. CONCLUSION AND DISCUSSION 

The matrix effective Hamiltonian and the matrix effective Lagrangian were constructed in the product integral. 
The matrix operator approach in the product integral formulation provided a unified picture to quantum evolution 
operators for parameter-dependent quantum systems. In the matrix effective Hamiltonian approach, the Schrodinger 
equation becomes the matrix equation (jl3p under the unitary transformation (jl2[) . The matrix equation can be 
solved to find the energy eigenvalues and the wave functions. As an application of the matrix effective Hamiltonian, a 
system of a heavy particle coupled to a light particle of harmonic oscillator (jisp was considered. The two-component 
wave function (j25p and the superadiabatic expansion (j26p were obtained in the product integral. The superadiabatic 
expansion can give asymptotic forms of the wave functions. 

On the other hand, in the matrix effective Lagrangian approach, the full time-evolution operator kernel (1361) and 
the matrix effective action (1381) were defined, and the Magnus expansion and the Kubo's cumulant expansion gave 
systematically the higher order corrections to the adiabatic kernel and the effective action. As a by-product, we 
found that up to the lowest Magnus expansion and the Kubo's cumulant expansion the effective action ([54| for the 
system (1401) of a heavy particle coupled to a light particle of harmonic oscillator involves the new (mode-dependent) 
gauge potential ([55]) that has been known not to exist for the nondegenerate case (real eigenfunctions). The real 
eigenfunctions ((TT)) were chosen, and the induced gauge potential (l53l) has its origin in a sequence of transition 
accompanying the creation of two quanta followed by the annihilation of two quanta or vice versa in the instantaneous 
number-state representation. 

The holonomy (15^ of the gauge potential (Berry's connection) could be evaluated using the Magnus expansion, 
and it was observed that the holonomy was related to the Schlesinger's formula for the gauge field tensor ((64)) . The 
Schlesinger's formula gave the exact geometric phase. The geometric transition amplitudes for quantum evolution 
operators were also formulated in the product integral and found to depend only on the contour integrals of a homotopy 
class around complex degenerate points. 

The matrix operator approach to quantum systems, which seemed especially applicable to parameter-dependent 
quantum systems, may provide an alternative method to perform further both the canonical quantum field theory 
and the path integral formalism. Moreover, the matrix approach is nonperturbative as far as the spectra resolution 
could be completely employed and the product integral of matrix operator could be found in a closed form. But the 
Magnus formal solution (expansion) for the product integral provides at least a perturbation method, which appears 
quite different from the conventional ones in that the first term in the Magnus expansion is the leading term as in Sec. 
mil The matrix effective Hamiltonian approach might be applied to canonical quantum gravity, whose wave function 
is a generalization of the adiabatic wave function to a time-dependent Klein-Gordon like equation. It may be used to 
find the complete set of solutions for such systems and provide the Hilbert space structure to the set of solutions [13] . 
The meaning of the new (mode-dependent) gauge potential for quantum gravity is left for a future investigation. 
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The matrix effective Lagrangian approach might be apphcable to finding systematically the back-reaction of subsys- 
tems. The higher order back-reactions are expected to contribute small amounts for the lower modes in the adiabatic 
regimes, in which the matrix effective Lagrangian method is efficient, but may contribute significant amounts for the 
higher modes. In the other case such as level crossings where the adiabatic method fails to work, the matrix-valued 
gauge potential contains a regular singular part and the matrix effective Lagrangian method in Sec. IIIII might not 
work. One needs to treat the singular part carefully and may expect geometric factors to enter in the path integral 
depending on the contours of a homotopy class of simple poles of the matrix-valued gauge potential as in Sec. I VII 
The singular gauge potential introduces geometric factor for the quantum evolution operator in Sec. IVIl Quantum 
theory of vortices in nonabelian gauge theory will provide such a challenging problem. 
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Appendix A: Magnus Formal Solution (Expansion) 



We briefiy summarize the Magnus formal solution (expansion) to the PI whose detailed description can be found 
in Appendix of Ref. [Tsj . An operator- valued first order linear differential (evolution) equation 



dt 

has a formal solution in the PI UM 



'""^^ = omit) (Ai) 



U{t) =Wciq>[0{t')dt']U{to). (A2) 



to 



The PI may be known in a closed form in a certain case; otherwise, we should rely on the formal (perturbative) 
solution, known as the Magnus expansion (l6| . The PI can be written in the Magnus expansion 

J|exp[0(t')rfi'] =exp[ / M[0{t')]dt']. (A3) 

to 

The commutator approach to the Magnus expansion follows from equating the derivative of the both sides 

M = ^-^(ad/ M[0{t')]dt'Yo{t), (A4) 

n=0 ^' •'^0 



where Bn are the Bernoulli numbers. Then the Magnus expansion can be found as an infinite sum [l 

»fl OO / t2 

- V . - V . ^ 

ni=0 \ •^*n ri2=0 \ „3=0 



M^Y.^nA^d Y.^nA^d (ad/ (•••) 0(^4) Oih)] 0{h)] Oih), (A5) 

„^^Q \ Jtn n \ Jtn n \ Jtn I I I I 

where 5* = Bn/n\ and 

(adA)"B= (adA)""^B], (adA^B^B. (A6) 

Here the integration are ordered integrations from the inside to the outside. One may insert an expansion parameter 
such as 1/h to the operator O and perform the multi-power expansions and the multi-partitions of integers into 
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nonnegative integers. The first few terms of the Magnus expansion are 

Af(o) = O, 

Af(3) = (-^)'(ad j (ad j (ad 1 o)o)o)o + (-1) l(ad j (ad 1 o)'o)o. (A7) 



The Magnus formal solution (expansion) \Kb\ can be used for the MEL in Sec. IIIII and the exponential form of the 
superadiabatic basis in Sec. IIVI since the commutators may get rid of large terms. 

There are other methods, such as the Burum's method z7] and the Salzman's method 28], which follow from 
the formal equivalence between the PI and the time-ordered integral. The Magnus expansion in these methods is 
an exponential operator in a series of combinations of time-ordered integrals. These methods might be used for 
calculating higher order geometric phase in Sec. IIVI 

We now list some useful properties used in this paper. The general perturbation formula (sum rule) is 

t t 
\{e^^[{Oi{t') + Oii{t'))dt'] ^P{t)\{ei,v[P'\t')Oii(t')P{t')dt'], (A8) 

tQ to 

where 

t 

P{t) = ^eyi^[Oi{t')dt'] . (A9) 

to 

The similarity rule for an invertible operator M{t) is 

J|exp[0(i')dt'] = M{t)\[e^^\(M-\t')0{t')M{t') - M-\t')^^^^ M-i(to). (AlO) 

to to 

The similarity rule resembles the diagonalization of the matrix (similarity transformation for a symmetric matrix) in 
a certain sense. 



Appendix B: Two-Component Wave Function 

We review and summarize the method to find the two-component wave function developed in Rcf. [TBj. The MEH 
equation ([M)) 

^-^^ = 2M[lVi{X) + Hii^oiX) ^ El]^iX), (Bl) 
can be written as a two-component equation of the first order 



d f ^{xn f A / *(X) 



dX \ ^ i " \2M{IVj{X) + Hji.d{X) -EI) Oj I 



The solution of the matrix differential equation (jB2|) may be given in the PI by 



, dX / X, 



J|exp 



/ 

2M{IVi{X') + HuAX') - EI) 



dX' 



\ dXi 



The exponent of the PI takes the form 



(B2) 



\YI^ e'A'ix')dx'y^ V l\ fU'' e^^'^^'^'^^' V . 

^Yi^ e'A'(x')dx'yij [2M{IVi + Hii,D-EI) 0) [ -q^ j^^^j 
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and with the aid of the similarity rule (jAlOl) . the two-component wave function can be written as 

„^A•(x')dx'^-l 



, dX , 



i^yx ^,A'(,x')dx'y^^ 



X 

n 







iA*{X') . , , 

2M{IVi + Hii^D - EI) iA*{X') ' 







Xo iA''(X')dX' 



n 

Finally, the two-component wave function is given by 



dXo , 



d^{x,X) 
dX 







<P'{x,X) 



X 



xn 

^mx',Xo)^ix',Xo)dx' 

m^^dx' 



I 



iA*{X') . . , 

2M{IVi + Hu,D - EI) iA*{X') ' ° 



/$(x',Xo) 



dXo 



(B5) 



(B6) 



Appendix C: Kubo Cumulant Expansion 

We summarize the first few terms of the cumulant expansion [25j 



{0)c = 


(0), 






{o\ = 










{0') 




■f 2(0)3, 




{0') 


-4(0)(03) 


- 3(0)2(02) + 12(0)3(0) - 6(0)4 
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